! just a test program

program untitled1()
   use VariableDef
   implicit none
   
   ! comment for ShockP, begin
   Tt=3000.0
   Pt=7000000.0
   g=1.2
   Rg=320.0
   p0=10000
   
   pi=3.141592653
   el=1e-5
   ep=1e-5
   erho=1e-5
   eV=1e-5
   etheta=1e-5
   delta=1
   icor=6
   
   yt=25   !0.025
   rtu=50   !0.050
   rtd=12.5   !0.0125
   thetaa=35*pi/180
   thetae=10*pi/180
   xe=250   !0.25
   ! comment for ShockP, end
   
   ! comment other than ShockP, begin
!   pi=3.141592653
!   ep=1e-5
!   delta=0
!   g=1.40
!   Rg=287.04
!   Vuu=1041.6
!   pu=100000.0
!   rhou=1.1613
!   Tu=300.0
!   Mu=3.0
   ! comment other than ShockP, end
   
   ! for ShockP, begin
!   x1=0.22113
!   y1=0.28317
!   V1=691.51
!   theta1=30.0*pi/180
!   p1=635590
!   rho1=3.6782
!   eps1=52.014*pi/180
!   x2=0.26991
!   y2=0.26988
!   V2=697.16
!   theta2=29.536*pi/180
!   p2=621320
!   rho2=3.6189
!   call ShockP()
   ! for ShockP, end
   
   ! for FreeP, begin
!   x2=0.34226
!   y2=0.12312
!   V2=2532.3
!   theta2=14.165*pi/180
!   p2=53164
!   rho2=0.12491
!   x3=0.32979
!   y3=0.12351
!   V3=2511.7
!   theta3=15.317*pi/180
!   p3=60000
!   rho3=0.13813
!   call FreeP()
   ! for FreeP, end
   
   ! for SymmetryP, begin
!   x1=0.079625
!   y1=0.001290
!   V1=2306.3
!   theta1=0.8866*pi/180
!   p1=170250
!   rho1=0.32947
!   x3=0.076195
!   y3=0.0
!   V3=2280.4
!   theta3=0.0*pi/180
!   p3=190770
!   rho3=0.36225
!   call SymmetryP()
   ! for SymmetryP, end
   
   ! for IWall, begin
!   x1=0.005495
!   y1=0.026020
!   V1=1728.9
!   theta1=24.090*pi/180
!   p1=1154500
!   rho1=1.6240
!   x3=0.005085
!   y3=0.026080
!   V3=1726.8
!   theta3=24.000*pi/180
!   p3=1160400
!   rho3=1.6309
!   x4=0.005283
!   y4=0.026170
!   theta4=25.000*pi/180
!   call IWall()
   ! for IWall, end
   
   ! for DWall, begin
!   x2=60.480
!   y2=59.625
!   V2=2274.2
!   theta2=30.122*pi/180
!   p2=195970
!   rho2=0.37046
!   x3=55.943
!   y3=58.845
!   V3=2252.9
!   theta3=30.752*pi/180
!   p3=214530
!   rho3=0.39947
!   call Wall2c()
!   call DWall()
   ! for DWall, end
   
   ! for InnerP, begin
   x1=2.7 !0.131460
   y1=0.01611 !0.040118
   V1=2218.1915 !2603.5
   theta1=0.0 !18.191*pi/180
   p1=240900 !34042
   rho1=0.44246 !0.086151
   x2=2.5877 !0.135683
   y2=0.02948 !0.037123
   V2=2178.1024 !2609.2
   theta2=0.008777 !16.422*pi/180
   p2=288251.1404 !32781
   rho2=0.5117 !0.083482
   call InnerP()
   ! for InnerP, end
   
   write(*,*) x4,y4,V4,theta4*180/pi,p4,rho4
   write(*,*) 'complete !'
end program untitled1
